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Abstract. It has been pointed out recently that the quadrupole-octupole alignment in the 
CMB data is significantly affected by the so-called kinetic Doppler quadrupole (DQ), which 
is the temperature quadrupole induced by our proper motion. Assuming our velocity is 
the dominant contribution to the CMB dipole we have vjc = /3 = (1.231 ± 0.003) x 10“^, 
which leads to a non-negligible DQ of 0(/?^). Here we stress that one should properly take 
into account that CMB data are usually not presented in true thermodynamic temperature, 
which induces a frequency dependent boost correction. The DQ must therefore be multiplied 
by a frequency-averaged factor, which we explicitly compute for several Planck CMB maps 
finding that it varies between 1.67 and 2.47. This is often neglected in the literature and 
turns out to cause a small but non-negligible difference in the significance levels of some 
quadrupole-related statistics. For instance the alignment significance in the SMICA 2013 
map goes from 2.3 cj to 3.3 cj with the frequency dependent DQ, instead of 2.9 cj ignoring the 
frequency dependence in the DQ. Moreover as a result of a proper DQ removal, the agreement 
across different map-making techniques is improved. 
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1 Introduction 

Since the observation of the anisotropies of the Cosmic Microwave Background (CMB) by 
COBE and WMAP satellites several authors have pointed out the existence of possible 
anomalies with respect to the standard paradigm of an isotropic ACDM cosmological model 
with Gaussian adiabatic and nearly scale-invariant primordial perturbations (see [1, 2] for an 
overview). The first task when investigating possible anomalies to any well-established model 
is to ensure that all systematics are taken into account. 

Among the many systematics involving a construction of a CMB map we focus here on 
those affected by our peculiar velocity with respect to the CMB, which induces both a Doppler 
and an aberration effect. The former’s main effect is to introduce a large dipole. But the 
effects of a boost to a CMB map extend well beyond the dipole: aberration changes the arrival 
directions of photons and Doppler introduces a multiplicative factor and as a consequence both 
effects correlate the spherical harmonic coefficients in a non-trivial way (see for instance [3- 
5]). It was in fact realized in [5-7] that Planck would be able to detect this signature and 
obtain an independent measurement of our velocity, which was subsequently done in [8]. The 
velocity found this way is in agreement with the one responsible for the dipole, reinforcing its 
interpretation as having a non-primordial origin. Moreover, our velocity influences at least 3 
of the discussed CMB anomalies: the low-quadrupole, the hemispherical asymmetry and the 
quadrupole-octupole alignment, all of which we now briefly discuss. 

The first one is part of the observed lack of power in the largest scales of the sky. It has 
been noted since COBE that assuming a power law spectrum the first multipoles of the CMB 
are rather low [9]. This remained true with WMAP [10-14] and Planck [15], and the spectrum 
for £ < 40 is currently in tension with the higher multipoles at around 2.5-3 it. The amplitude 
and spectral index are in fact basically fixed by fitting the high multipoles while the \ow-^ 
are consistently below the average value. Among these large scales one of the most debated 
one is the quadrupole [11, 16, 17] whose amplitude is ~ 200(^iP)^, almost 6 times lower than 
the theoretical expectation for Planck’s ACDM parameters, which is around 1170(;uiP)^ [18]. 
Note however that the Doppler effect, in addition to a large contribution to the dipole, also 
produces a small kinetic contribution to the quadrupole, which has to be taken into account 
properly. Moreover, it was realized in [19] and later also in [20, 21] that such a quadrupole 
contribution has a frequency dependence, which makes its subtraction more subtle. We will 
address this issue quantitatively later on in this short note. 

The hemispherical asymmetry is related to an apparent modulation of the power of the 
CMB anisotropies in the sky. In particular, by splitting the CMB in two hemispheres there 
seems to be in WMAP [22-26] and also in Planck [2] an indication for a power asymmetry 
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along the direction {l,b) = (225°,!°). In vl of [2] such asymmetry was claimed to extend 
with a large amplitude up to £ = 1500; at the same time, it was also noted in [8] that the 
Doppler effect due to our motion should induce a 0.25% power asymmetry. Later we have 
shown quantitatively in [27] that both Doppler and aberration induce a power asymmetry 
and that it grows at small scales, reaching 0.6% at high i and that the subtraction of such 
effect is crucial in the estimation of the asymmetry if considering I up to about a thousand. 
In fact in revisions of [2] and in [28] it was shown that the significance at high I disappears 
after such subtraction so that a scale-independent power asymmetry is now actually excluded 
by Planck and a dipolar modulation of power can survive at most up to £ < 600, though 
with a marginal significance of less than Su (which also agrees with earlier results in [29] for 
WMAP). 

Finally, it seems that the three lowest multipoles might possibly not be randomly 
distributed in the sky as there seems to be some alignment between them in the WMAP 
data [12, 30]. In fact a direction for the quadrupole and octupole can be defined either by 
maximizing the analogous quantity to the angular momentum [31] or by using the so-called 
multipole vectors [32]. In both cases they are found to be extremely close to each other, which 
is unlikely in a Gaussian isotropic distribution: the significance of such an alignment has been 
estimated to correspond to a p-value of about 0.1% (see e.g. [33]). With the release of the 
Planck data such an alignment seemed at face value to have turned to be less significant [2] 
to a p-value of 1% — 2%, but it has been pointed out [33] that this has to be corrected for the 
Doppler effect, which induces a velocity dependent quadrupole (kinetic Doppler quadrupole, 
DQ) and this brings up again the significance of the alignment to a level similar to the one 
found in WMAP of about 0.05% — 0.3%; note moreover that such common direction seems 
to be also correlated with the dipole direction, giving similar significance levels for a statistic 
which involves 1" = 1, 2, 3 (see Table 5 and 7 of [33]). 

However, we stress here that due to the way in which the CMB data are presented, both 
by WMAP and Planck, such a correction is actually frequency-dependent, as it has been 
already shown in [19] and similarly to what has been discussed in [8] in a different context. 
This results in an effective multiplicative correction factor for each frequency map. This factor 
is very different from 1 for the single frequencies, being especially large at high frequencies, 
and therefore it has a potentially large effect in a linearly combined temperature map, such 
as the Planck SMICA and NILC maps. To our knowledge such a frequency dependent factor 
for the DQ has only been quantitatively considered in [2], and it has been estimated to be of 
1.7 for both the SMICA and NILC maps and then it has been applied to the Wiener filtered 
maps of the Planck collaboration to show that the alignments have a significance of about 
1%. Other groups [33-36] have used instead just a factor of 1, which is significantly off (in 
one map the real value it is over 2.5, see below) and leads to a non-negligible modification of 
their results. In this note we make an independent estimate of this factor for several maps, 
including two which are not constructed by the Planck collaboration. 

2 Frequency dependence of DQ in CMB Experiments 

Both WMAP and Planck collaborations are showing at each frequency intensity maps which 
are assumed to be proportional to the CMB temperature, through a linear series expansion. 
However the exact relation between the intensity measured in our frame and the rest-frame 
CMB temperature contains a frequency dependent factor. Let us briefly recall how such a 
factor arises for the DQ, similarly to the treatment of [19]. 
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A boost with a velocity 'v/c = [3 has two effects on an image of the sky: Doppler effect 
and aberration. The first effect is a multiplicative direction dependent factor while the second 
effect changes the apparent arrival direction of photons [37]: 


r'(n') = 


r(n) 


7(1 - /3 • n ) 


( 2 . 1 ) 


where n' is the (aberrated) arrival direction of the photons and r'(n') is the temperature 
observed in our reference frame. Here we assume the CMB temperature in the rest frame 
T(n') to be the sum of a homogenous component plus perturbations in this way: 


r(n) =ro + e<5r(n), 


( 2 . 2 ) 


where by convention we assume e = 10“® and so 5T{fi) is of order 1. 

If a map is completely homogeneous (e = 0 ) aberration has no impact, and so its leading 
effect is of order e (3. However the Doppler effect is non-trivial even if a map is completely 
homogeneous and in fact it induces a dipole of order (3 and an n**^-pole of order /?"■, as can be 
seen by expanding the multiplicative factor in eq. (2.1). Therefore, if our velocity is /3 ~ 10“^, 
there is an induced dipole of order 10“^ and a kinetic Doppler quadrupole (DQ) correction 
of order 10“®, which is not negligible. All other effects are at most of order of /3e or and 
therefore negligible in the discussion of the alignments (but important for other purposes, 
such as measuring our velocity through non diagonal correlations between multipoles, as it 
has been shown in [5-8]). 

As we already mentioned, WMAP and Planck frequency maps [38-41] are not produced 
in terms of temperature, but in terms of intensities multiplied by a proportionality factor and 
so the boost induces a frequency-dependent correction. An observer in our frame sees in fact 
an intensity at a given frequency u' given by: 

JV) = Jw(-)"= if” ■ (2.3) 

\i^ J eT(ft) _ 1 


where = z^^ 7 (l — /3 • h'). In the /3 = 0 limit the two frames coincide and fluctuations in 
Intensity are given at first order in e by: 


61 {u) 


2z/'^e^o 


T^ie^o -1 


ye5T{n) = Ke 


6T{h) 

To 


(2.4) 


The above approximate equation is used by the WMAP and Planck collaborations to 
define temperature as 5I{i')/K, which differs from the real thermodynamic temperature. We 
will refer to 6I{v)/K here as linearized temperature. We can in fact now write down the 
expansion at second order in f3 and first order in e of the linearized temperature: 


K 


6T(h) 

£ — 

To 


+ fiz + (3^z^Q{v') 



(2.5) 


where we have discarded terms of order (3 e 01 higher and where z = f3 ■ n and 


QW) = coth ( -— 


2T, 


2To 
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in agreement with [19].^ So, in addition to a dipole correction we also have a freqnency 
dependent (qnadrnpole) correction and constant shift to the monopole. Using the fact 
that To = 56.8 GHz we get that the correction factor Q(i^) has for Planck a valne of abont 
1.5 for the 143 GHz channel and of abont 2 for the 217 GHz channel, which are the main 
two channels which contain cosmological signal. However, in order to remove foregronnds 
all CMB maps come from combinations of more freqnencies, so we perform here an explicit 
compntation of the overall effective Qeff for specific maps. 

In this work we analyze 4 different CMB maps: 2 prodnced by the Planck Collaboration 
(SMICA 2013 and NILC 2013) and 2 prodnced by an independent gronp [42]. For both 
SMICA 2013 and NILC 2013 Planck provided inpainted maps where the inpainted regions 
are generated with a constrained Ganssian realization which assnmes isotropy [43]. Althongh 
it consists of only a small 3% fraction of the sky, it is interesting to assess the effect of the 
inpainting assnmptions, especially since the fidncial cosmology assnmed in the inpainting is 
not discnssed in detail in [39]. Using the pnblic Planck inpainted maps [44] we find resnlts 
that agree with the ones in [33], where a different inpainting techniqne was nsed. 

It shonld be noted that the inpainting masks are smaller than the so-called confidence 
masks of Planck, which cover 11% and 8% of the sky for SMICA and NILC, respectively. So 
there conld be some residnal foregronnd contamination on these inpainted maps. Probably 
for this reason, Planck collaboration also analyzed an alternative SMICA map with a mnch 
larger 27% mask, inpainted with a Wiener filter techniqne [2]. We conld not compnte the 
significance of this map becanse it is not pnblicly available. Instead, for SMICA and NILC 
we consider also the alternative maps provided in [35] , the only two differences being that the 
masked region are inpainted with a different techniqne and and that for NILC the inpainted 
region correspond to the 8% confidence mask instead of the 3% one. We dnb these maps 
SMICA (R14) and NILC (R14). 

In what concerns the qnadrnpole and octnpole, all 4 maps can be snmmarized as a simple 
linear combination of the different freqnency channels with a set of weights, which snm to 
1. In all cases except one the weight-vectors contain 9 entries and make nse of all 9 Planck 
channels; the remaining map has a weight-vector containing 14 entries, dne to inclnsion also 
of WMAP’s 5 channels. Once we have the weights , we jnst compnte: 

= 'Y^Q{yi)wv,- (2.7) 

i 

We will qnote the weights in order of increasing channel freqnency, so that the 5 WMAP 
entries correspond to the freqnencies (23, 30,40, 60, 90) GHz and the 9 Planck ones to the 
(30,44, 70,100,143, 217,353, 545, 857) GHz channels. 

Given these weights, there is an important caveat: the low freqnency channels (the ones 
measnred by the LFI instrnment) of the 2013 Planck release have been already treated in 
snch a way to remove the DQ, bnt a bng in the code removed only half of the freqnency- 
independent effect (see the Appendix of [45]). Therefore for the LFI channels we apply 
Q{i') —)• Qii') — 0.5. The high freqnency channels (HFI) instead have not been treated at all, 
as described in Appendix of [46]. 

For WMAP alone the varions channels are all at qnite low freqnency and so the Q factors 
are between 1.01 and 1.22. In this case the snbtraction with Q = 1 obtained in [33] gives 
already a good approximation of the correct resnlt. 

^Note that there is a typo in the derivation of Eq. (2) of [19]: = "fv'{l + j3fi) should read v = 'yi/fl + Pjj,'). 
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For the Planck SMICA map, the weights provided in both Figures D.l of [39, 41] are 
not the hnal numbers we are interested in, because they are in fact a product of 3 different 
functions: the actual CMB temperature weights, the conversion factors from Antenna Tem¬ 
perature to linearized temperature (he. the ratio between K of eq. (2.4) and its Rayleigh-Jeans 
small u limit) and a beam and HEALPix^ [47] pixel window function correction [39]. For the 
quadrupole and octupole, the last term is irrelevant. We obtained the linearized temperature 
weights by multiplying the unit conversion factors to the carefully inspected values of Figure 
D.l of [39] 

WSMICA 13 = {0.01, -0.094, -0.113, -0.256,1.184,0.436, -0.174,0.005, -8 x 10"®} . (2.8) 

As a consistency check such numbers now sum to 1 to good precision, as they should. Keeping 
in mind the different units, wsmicais also agrees very well with the values read independently 
in a similar way by [48]. We could do the same with SMICA 2015, but it is unclear how the DQ 
was treated there. In fact in [41] it is stated that the DQ has already been subtracted and this 
is clearly visible in the difference of the 2013 and 2015 maps, but it is not clearly specihed 
if the correct Q(u) factors have been applied to the single frequency channels. Moreover, 
information regarding the calibration of the 2015 LFI channels has not yet been released. 

For the Planck NILC map only the hrst needlet band, which has a constant weight across 
the sky, is relevant. The weight vector for NILC 2015 map was obtained by careful inspection 
of Figure B.2 in [41] obtaining: 

WNILC15 = {0.05, -0.055, -0.7,0.25,1.51,0.03, -0.09, 0.005,0} . (2.9) 

Unfortunately the weights for NILC 2013 are not reported in the Planck papers, so we cannot 
conhrm the value of 1.7 obtained in [2]. In fact, the value obtained in 2015 differs from 1.7, 
which indicates that the weights changed between both releases. Thus, for the particular case 
of NILC we use this quoted value from Planck. 

The LGMCA maps [42] uses a component separation method which is claimed to be 
able to reconstruct the full sky without the need of a mask nor inpainting. Specihcally the 
authors have released two maps: the PRl map, which uses the 9 Planck channels and the 
WPRl which uses both the WMAP 5 channels and the Planck 9 channels. The corresponding 
weights could be obtained from the map-construction code which is publicly available [49]: 

wpRi = {0.087,0.087, -1.186, -0.668, 2.093, 0.885, -0.306, 0.008, -7 x 10"®} , (2.10) 

and 

wwPRi.WMAP = {0.062,-0.277,0.227,0.191,0.413}, (2 11) 

WWPRI,Planck = {0.147,-0.45,-0.827,0.254,1.412,-0.064,-0.09,0.002,6 x 10“®}. 

Note, again, that such PRl and WPRl weights respectively sum to unity. In such case we 
hnd that the factor is 2.47 and 2.05, respectively for PRl and WPRl. 

The hnal correction factors at i = 2 range from 1.67 to 2.47 and are given in Table 1, 
together with the change on the signihcance of the alignments, using several statistics. It is 
important to stress that individual Q(u) for the higher frequencies have much larger values 
for Planck: for instance they are 3.1 and 4.8 for the 353 GHz and 545 GHz Planck channels, 
and Qeff is the result of cancellations among large contributions. Note also that we have 

^http://healpix.sourceforge.net/ 
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checked that in the new release of Planck 2015 the weights have changed by a large amount, 
but as we already said we do not know if the correct overall Qefr has been already applied or 
not to the maps; this will probably be clear once the LFI calibration paper is released. 

Note that in principle Q{y) should be averaged over the bandpasses for each channel 
with an appropriate weight. As a simple check we have performed a weighted average with a 
step function between some minimal and maximal frequency. We have found values for the 
minimal and maximal frequency for each HFI channel in Table 3 of [50] and for LFI in [51] 
(although for LFI this is completely irrelevant). We find a 1% difference in the final Qeff) with 
respect to our previous procedure, so we have neglected such corrections in what follows. 


3 Summary of our Results 


Once we have the correct Qgff factors we have then computed the significance of the alignments 
using three types of statistics, which make use of the multipole vectors (see [33] and references 
therein). The multipole vectors are a different way of expressing the same information as the 
a^m, which consists in one amplitude and I unit vectors (which sums to 2.^+1 real quantities). 
From the multipole vectors one can build the area vectors Wj, which are cross products of the 
multipole vectors. One can also define instead directions fii of maximal “angular momentum” 

^^771 |ci7m| , (3'1) 

m 


one for each i [31]. For £ = 2 in fact n 2 = W 2 /IW 2 I, but for higher multipoles the relationship 
is more complicated; nevertheless they both probe similar effects [52]. From the computational 
point of view, however, it is much simpler to compute the multipole vectors, as these can be 
found simply by computing the roots of a polynomial built using the a^m’s of a given i [53] . For 
hi, on the other hand, the standard practice is to rotate the in many directions (often 
millions) on the sky and compute explicitly the “angular momentum” in each case. In order 
to compute h^ we made use of the code developed in [35] and made publicly available [49]. 
We also cross-checked that n 2 = W 2 /IW 2 I was found to hold comparing both techniques. 

The first statistics is the simple relative angle between n 2 and h^, where we assign a 
probability assuming a flat distribution of n 2 ■ h^, which corresponds to random orientation 
of the two unit vectors. The other two statistics that we used are the S and T variables [33] 
which are defined instead with respect to a particular direction h and they involve the area 
vectors of the quadrupole and of the octupole as follows: 


= 

n 


i=i 


T = 1 - - 
n 


n 


-A^ 


J=1 


(3.2) 


where Aj = {wj ■ h\. In the present treatment we will only use as a specific direction n the 
direction of the dipole, which has been shown to be the one with most significant alignment 
according to the analysis of [33]. In other words we test the hypothesis that the dipole, 
quadrupole and octupole are all aligned. 

In Figures 1 and 2 we depict the variation of the quadrupolar angular momentum as a 
function of the direction used to expand the spherical harmonics. The direction that maxi¬ 
mizes this value defines n 2 . In both figures, we show on the left hand side the results ignoring 
the DQ and in the right the properly DQ corrected results. Note that the right-hand plots 
are in better agreement and that in the PRl case, which is most affected by the kinetic 
quadrupole (Qeff = 2.47), n 2 changes substiantially. 


6 



Figure 1. Total quadrupolar angular momentum - defined by Eq. (3.1) - for different directions 
of the spherical harmonic expansion, for the SMICA 2013 map and galactic coordinates. Red (blue) 
regions denote the highest (lowest) values. The red-most point correspond to the n 2 direction. The 
actual angular momentum values (i.e. the scale of the plot) is not relevant here. Left: raw map 
(Qeff = 0). Right: map corrected for the DQ (Qeff = 1.67). 



Figure 2. Same as Figure 1 for the LGMCA-PRl map (and Qeff = 2.47). Note that the difference 
in n 2 in this case is much larger, although the overall change in the distribution is not very big. 


We will generally quote our results in a levels (instead of p-values) which are defined 
through the Gaussian distribution: 

na = -v/2Erf“^(l — p-value) . (3-3) 

In other words, we use just a straightforward generalization of the commonplace association 
of 1, 2, 3cr as shorthand notation for 68.3%, 95.4%, 99.73%, (independently of whether the 
underlying distribution is really Gaussian) and so forth. 

The results for our three alignment statistics are listed in Table 1. We compare the values 
in 3 cases: without correcting for the kinetic quadrupole; correcting with the wrong Qeff = Ij 
as done in [33-36]; and finally with the correct result, using the appropriate Qeff for each map. 
These results are based on a 1-tailed statistical test. In Table 2 we list the corresponding 
two-tailed statistics for the case of the appropriate QeS- In fact, if one considers a quadrupole 
almost exactly orthogonal to the octupole also as an anomaly, then one must use both tails of 
the distributions for n 2 -n^. A similar concern applies to S and T, whose distributions (taken 
from interpolation of Figure 6 of [33], realized with 10® inpainted simulations) are peaked 
around a central value in an almost Gaussian way and so both tails should be considered to 
be anomalous.^ 

^It is maybe instructive to pose the following sociological question: had the dipole, quadrupole and octupole 
been observed to be almost exactly perpendicular to each other, would there have been many papers in the 
literature studying this phenomena? If one thinks this would have happened, the 2-tailed distribution is the 
correct assessment of the significance as far as the Cosmology community is concerned. 
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Map Qeff ^2 • 1^3 (prob) S (prob) T (prob) 


Uncorrected values 


SMICA 2013 

0 

0.9811 

(2.3cj) 

0.742 

(2.8a) 

0.923 

(2.9a) 

SMICA 2013 (R14) 

0 

0.9833 

(2.4a) 

0.744 

(2.8a) 

0.926 

(2.9a) 

NILC 2013 

0 

0.9901 

(2.6cj) 

0.734 

(2.7a) 

0.922 

(2.9a) 

NILC 2013 (R14) 

0 

0.9853 

(2.4a) 

0.736 

(2.7a) 

0.925 

(2.9a) 

LGMCA-PRl (R14) 

0 

0.1739 

(0.2a) 

0.558 

(1.4a) 

0.691 

(1.0a) 

LGMCA-WPRl (R14) 

0 

0.8851 

(1.6a) 

0.696 

(2.4a) 

0.873 

(2.3a) 


Boost-corrected values, incorrect QeS 


SMICA 2013 

1 

0.9964 

(2.9a) 

0.775 

(3.1a) 

0.937 

(3.1a) 

SMICA 2013 (R14) 

1 

0.9973 

(3.0a) 

0.778 

(3.1a) 

0.940 

(3.2a) 

NILC 2013 

1 

0.9987 

(3.2a) 

0.777 

(3.1a) 

0.944 

(3.3a) 

NILC 2013 (R14) 

1 

0.9997 

(3.6a) 

0.772 

(3.0a) 

0.943 

(3.3a) 

LGMCA-PRl (R14) 

1 

0.9824 

(2.4a) 

0.648 

(2.0a) 

0.838 

(2.0a) 

LGMCA-WPRl (R14) 

1 

0.9792 

(2.3a) 

0.747 

(2.8a) 

0.906 

(2.6a) 


Boost-corrected values, correct Qeff 


SMICA 2013 

1.67 

0.9990 

(3.3a) 

0.788 

(3.2a) 

0.940 

(3.3a) 

SMICA 2013 (R14) 

1.67 

0.9985 

(3.2a) 

0.791 

(3.2a) 

0.943 

(3.3a) 

NILC 2013 

1.7 (*) 

0.9974 

(3.0a) 

0.794 

(3.2a) 

0.949 

(3.4a) 

NILC 2013 (R14) 

1.7 (*) 

0.9990 

(3.3a) 

0.786 

(3.2a) 

0.947 

(3.3a) 

LGMCA-PRl (R14) 

2.47 

0.9844 

(2.4a) 

0.752 

(2.9a) 

0.927 

(2.9a) 

LGMCA-WPRl (R14) 

2.05 

0.9975 

(3.0a) 

0.806 

(3.4a) 

0.948 

(3.4a) 


Table 1. 1-tailed statistics for the quadrupole-octupole alignments. For SMICA we use both the 
public map and the one in [35], to depict the effect of different inpainting techniques. The use of a 
single tail of the distributions is the correct approach if one considers a quadrupole almost exactly 
orthogonal to the octupole (and also to the dipole, in the case of S and T) as not anomalous. For 
QeS = 1 we reproduce the results in [33, 36[. (*) For NILC 2013 we were unable to compute QeS as 
the different needlet weights were not made available. We instead rely on the number quoted in [2[. 


Having computed the relevant Qeff also has an impact on the amplitude of the quadrupole 
and so we also show results for the “low quadrupole” anomaly. We compute the significance 
assuming as usual a distribution with 5 degrees of freedom and centered around the 
theoretical estimate, computed using the best-fit ACDM parameters. For Planck SMICA 
2013, this amounts to 1170 We list the results in Table 3, both using a 1-tailed 

and 2-tailed statistics. It is clear that the quadrupole is in itself only at most marginally 
anomalous; in the case of allowing also the possibility of having a high quadrupole (which we 
see no reason not to) it is not anomalous at all. The effect of DQ correction turns out to be 
an overall increase of around 0.2a in all the significance levels in both statistics. This seems 
to be of the same order as the difference among the different map making techniques, which 
differ mostly in their treatment of the regions close to the galactic plane. 

We now also briefly discuss the possibility of other astrophysical (non-primordial) effects 
contributing to the quadrupole and octupole. Two such effects were studied in [35]: the ISW 
contribution (which affects the low-£ modes of the CMB) and the kinetic Sunyaev-Zel’dovich 


- 8 - 









Map 

h 2 ■ ns (prob) 

S (p 

>rob) 

T (prob) 

Boost- 

corrected values, 

correct 

Qeff 



SMICA 2013 

0.9990 (3.1ct) 

0.788 

(3.O0-) 

0.940 

(2.6(t) 

SMICA 2013 (R14) 

0.9985 (3.0ct) 

0.791 

(3.O0-) 

0.943 

(2.7(t) 

NILC 2013 

0.9974 (2.80-) 

0.794 

(3.1a) 

0.949 

(2.80-) 

NILC 2013 (R14) 

0.9990 (3.10-) 

0.786 

(3.O0-) 

0.947 

(2.80-) 

LGMCA-PRl (R14) 

0.9844 (2.20-) 

0.752 

(2.60-) 

0.927 

(2.4a) 

LGMCA-WPRl (R14) 

0.9975 (2.80-) 

0.806 

(3.20-) 

0.948 

(2.8a) 


Table 2. 2-tailed statistics for the quadrupole-octupole alignments, using the appropriate The 

use of both tails of the distributions is the correct approach if one also considers a quadrupole almost 
exactly orthogonal to the octupole. In the case of S and T the distribution is bell shaped and so both 
small and large values are to be considered anomalous. 


Map 

raw C 2 (piAT)^ 

corrected C 2 

1-tail prob 

2-tail prob 

SMICA 2013 

251 

228 

2.1a 

1.1a 

SMICA 2013 (R14) 

235 

214 

2.2a 

1.2a 

NILC 2013 

219 

196 

2.2a 

1.2a 

NILC 2013 (R14) 

230 

209 

2.2a 

1.2a 

LGMCA-PRl (R14) 

284 

224 

2.1a 

1.1a 

LGMCA-WPRl (R14) 

189 

158 

2.4a 

1.4a 


Table 3. 1 and 2-tailed statistics for the low-quadrupole, using the appropriate The use of 

both tails of the distributions is the correct approach if one also considers a high-quadrupole as an 
anomaly. The significance of the low-quadrupole in itself is marginally significant only if one neglects 
this. 


effect (kSZ) from CMB photons scattering off moving free electrons in the local group. They 
found that the former made a significant impact on these large scale-modes, whereas the latter 
effect was too small. For the low-quadrupole signihcance the net effect was small because ISW 
subtraction affects both the measured value and the theoretical expectations. For the i = 2, 3 
alignments, however, they were claimed to become insignihcant after ISW corrections. This 
would imply that the alignments, if not a statistical fluke, might be fundamentally due to 
some interactions of the photons with large-scale structure instead of due to some primordial 
physics. Nevertheless the amplitude of the ISW effect was estimated in [54] to be signihcantly 
smaller and in agreement with the theoretical estimate of [55], which if conhrmed would make 
its effect on the alignments smaller than the estimate of [35]. 

4 Conclusions 

We have analyzed the impact of the Doppler effect on the quadrupole for WMAP and Planck 
and we have found that it induces a correction factor which is frequency dependent and large 
at high frequencies. While for WMAP, which operated at smaller frequencies, such a factor 
is of small importance, for Planck this is not the case. 

The effective correction factor for Planck is equal to: 1.67 for the SMICA 2013 map, 1.7 
for NILC 2013, and 2.47 and 2.05 for the PRl and WPRl full-sky reconstructions with the 
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LGMCA component separation method of [35]. We have applied such a factor to the analyses 
of [33, 35] hnding a variable change in the signihcance of, in some cases as large as Icr. 

After properly correcting for the DQ we hnd that the quadrupole-octupole alignment 
in the released inpainted versions of SMICA 2013 and NILC 2013 maps (which replaces 3% 
of the sky) as well as the full-sky WPRl map are signihcant enough to be taken in serious 
consideration as a possible challenge to ACDM. These results are also found in alternative 
inpaintings of SMICA and NILC, the latter of which masks a larger 8% portion of the sky. On 
the other hand the full sky PRl map and the values quoted in [2] (based on different versions 
of SMICA and NILC, which are not public, treated with an inpainting performed on much 
larger region of the sky, to wit 27%) yield lower values, even after the correct subtraction of 
the DQ. They lead to 2.2cj, 2.4a and 2.3a (2-tailed statistics), respectively (where we made 
direct use of the values quoted in [2]). This emphasizes the need for a better understanding 
of the assumptions and procedures of how the galactic region is treated, i.e. by masking and 
inpainting with various assumptions or by direct reconstruction. 

The frequency dependency of the quadrupole has been so far mostly neglected in the 
literature. We have shown in fact in [56] that it was also ignored in the calibration procedure 
of Planck 2015 data (which is based on the orbital dipole) even though it is of the same order 
of magnitude as other effects considered. An inclusion of the Q(i^) term could thus improve 
the calibration, for instance suppressing the leakage of temperature into polarization modes. 

It is important to realize that the CMB largest scales are completely cosmic-variance 
dominated and so, decreasing the experimental noise yields no improvements. These can 
come only through better understanding of the systematics, reduction of the masked portion 
of the sky or measurements of large-scale polarization. In effect, the DQ is also relevant for 
polarization, so its proper treatment will also be important there. 

It is also relevant to stress that even a small change in the weights of the different 
channels would result in a potentially large change in the factor Qgg and therefore in the 
direction of the quadrupole. So, in the presence of signihcant uncertainties in foreground 
subtraction, which would change the relative weights, one has to be very careful in applying 
the correct Qeff factor to the quadrupole. 
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